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ABSTRACT 

This thesis reports the results of an experimental investigation into 
the spanwise correlation of the fluctuating longitudinal velocity component 
in the early wake of a flat plate placed parallel to a uniform stream and 
allowed to vibrate about a vertical axis through the leading edge. The 
experimental program provides for the measurement of fluctuating total head 
signals at one fixed and one variable spanwise position on the wake center- 
line one body thickness downstream of the trailing edge. Analyzation of 
these signals gives the spanwise correlation of the wake structure expressed 
as a complex number and the power spectral density of the longitudinal 
velocity fluctuation along the span. In addition the mean local wake 
velocity is measured for each value of the ambient flow velocity. 

Ambient velocities of 8, 11, and 14 fps and a single plate of £" thick- 
ness and 2" chord, were used for all the tests which were carried out in the 
7 ^tt x q»» test section of an open-circuit water tunnel in the Hydrodynamics 
Laboratory at MIT. 

An equation is developed which relates the vortex-induced forcing 
moment on the foil to the measured structure of the early wake. Using this 
model, forcing moment coefficients are evaluated as a function of a "self 
excitation" parameter consisting of the ratio of the mean transverse velocity 
of the trailing edge to the free stream velocity. The role of spanwise 
coherence of the wake is examined with respect to its effect on the coupling 
between body vibration and forcing moment. 

Major Conclusion of this study is : 

The increase in the forcing moment coefficient, m^ , accompanying an 

o 

increase in the self excitation parameter is due to increasing correlation 
or coherency of the wake along the span rather than to increasing strength 
of the individual vortices being shed. 
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I INTRODUCTION 



A. Flow-Induced Vibrations 

Structures subjected to a flow field may, under certain circumstances 
experience flow- induced vibrations. These vibrations are undesirable and 
indeed may cause complete collapse of the structure due to fatigue induced 
by the vibrations. 

Unstreamlined or '’bluff" bodies in a flow field shed discrete vortice 
in the wake. The shedding of these vortices gives rise to periodic lift 
forces in a direction transverse to the flow direction. It has been shown 
in reference (2) that for bodies of finite thickness this vortex- induced 
lift force increases linearly with the transverse vibrational velocity, 
thus producing what may be called a "self-excited" vibration. The fluid 
mechanism involved in this feedback between motion and forcing force is 
unknown but may be expected to be due to an increase in the strength of 
the discrete vortices being shed and/or to an increase in the spanwise 
correlation of the vortices being shed. 

B. Scope of the Investigation 

The purpose of this investigation is to determine the cause of the 
amplitude dependence of the forcing moment so that engineers may devise 
steps to prevent the buildup of large amplitude vibrations. 
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II EXPERIMENTAL EQUIPMENT 

A. Water Tunnel 

The experiment reported herein was carried out in the closed-jet, open 
circuit water tunnel of the Hydrodynamics Laboratory at MIT. The water 
tunnel is fully described in reference (1). Figure 1 shows the water tunnel 
test section. Pump number III (rated capacity 3000 G.P.M. at a head of 45 
ft. and 1175 R.P.M. ) was used. 

Flow measurement was obtained using the tunnel contraction as a flow 
meter. Its calibration was checked and found to be the same as reported in 
reference (2). 

B. Test Plate 

The test plate used was plate D from reference (2), having the 
geometry shown in Figure 15* The plate was mounted vertically in the test 
section in such a way as to permit rotation about the center of its semi- 
circular leading edge. The plate is fabricated of aluminum. 

C. Plate Motion Measurement 

1. Strain Gages 

Tatnall type C9-121-R2C strain gages were cemented to the torsion 
spring which restrains the plate motion and were waterproofed. The four 
gages form a Wheatstone bridge which may be balanced by means of an external 
variable resistance. The bridge output voltage is proportional to the angle 
of plate rotation. See Appendix C of reference (2) for theory of calibra- 
tion of the strain gage bridge. 

2. Amplifier 

A decade amplifier model 220-A made by Ballantine Laboratories Inc. was 
used to amplify the strain gage bridge output. 
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3* Oscilloscope 

An Eico model 46 O DC-wide band oscilloscope was used for visual 
observation of the strain gage bridge output. 

D. Total Head Fluctuation Measurements 

1. Total Head Turbulence Gages 

The gages for measuring the total head fluctuations (turbulence) are 
fully described in references (3 and 4)« Figures 3 and 4, Appendix A, are 
photographs of the total head gages and Figure 15, Appendix A, is a drawing 
of the probes and test plate. 

2. Electrometer Preamplifier 

Two electrometers (one for each turbulence probe) were used to match the 
probe and amplifier impedences. 

The electrometers are fully described in references (3 and 4)» 

Figure 16 is a block diagram showing all electronic components. 

3. Amplifiers 

Two Dana model 2000 solid state, wideband, differential, low-level D.C. 
amplifiers with integral power supplies were used to amplify the turbulence 
probe signals. The amplifier has a flat response to ^3db from DC to 
approximately 70 KC. 

4. Oscilloscope 

A Dumont type 333 dual beam oscilloscope was used for visual observa- 
tions of the two turbulence probe signals. 

5« Portable Tape Recorder 

A Precision Instrument Co. portable tape recorder model PS 200 A was 
used to record the turbulence probe signals. 
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6. Endless Loop Tape Recorder 

A Sanborn model 2000 endless loop tape recorder was used to give a 
continuous input of the turbulence probe signals into the wave analysis 
system. 

7. Wave Analysis System 

To analyze the turbulence probe signals, a Technical Products Co. 

TP 625 Wave Analyzer System was used. This system consists of the following: 

a. One TP 626 Oscillator 

This oscillator produces an output frequency which is continuously 
variable from 97 KC to 122 KC over three different frequency bands: 

97.0 KC - 97.25 KC; 97.0 KC - 99-5 KC; and 97.0 KC - 122.0 KC. The 
oscillator output is used to modulate the signal to be analyzed in 
the TP 627 analyzers. 

b. Two TP 627 Analyzers 

These analyzers determine the frequency and amplitude of the 
individual components of a complex wave within the frequency range 
of 2-25»000 cycles. The results are measured in percent of total 
signal. 

These analyzers are equipped with two matched TP-218 type 
filters as required for operation in conjunction with a TP-645 
Multiplier. 

c. One TP-645 Multiplier 

This multiplier, when used in conjunction with the other components 
of the system, gives the cross power spectrum of two complex waves with- 
in the frequency range of 2-2 5 , 000 cycles. 
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III THEORY 

A. Forcing Moment Coefficient 

In reference (2) Noutsopoulos develops the equation of motion of a 
single degree of freedom plate spring system vibrating about an axis through 
the leading edge of the plate. The motion is caused by a forcing moment 
resulting from flow separation at the trailing edge of the plate. Assuming 
two-dimensional motion, the equation of motion of the plate may be written 
I t « + C of + K t ot = M 
where 1^ = the total moment of inertia 

= the total damping coefficient 
= the total stiffness of the system 
and c< = the angular plate displacement in radians cx - C 

From purely dimensional considerations, the simplest form of the 
forcing moment can be written 

M = m f \ p H(2b) z V 2 

where m^ is the spanwise average forcing moment coefficient, H is the span 

and V the ambient fluid velocity. When the forcing moment is periodic as 

in the case of excitation by vortex shedding, the moment coefficient is 

conveniently written in complex form as 

i(wt + 6 ) 
m f = m f e T ' 

o 

where w = angular frequency of vortex shedding from one side of the plate 
in cps 

= phase angle tetween the spanwise average forcing moment and the 
motion o<\ 



t = time in seconds 
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The moment on the plate may also be written in terms of a lift force 
and an appropriate moment arm where the lift force is a function of the 
circulation which is a measure of the vortex strength. 

Therefore 

M ~ 2bN 

where 2b is the chord length and 
N is the lift force . 



The lift force may be written as 



N = p V p H 

where T is the circulation on the plate and is a function of position along 
the span and time such that: 

f H 

N = p V I T (z,t)dz . 

0 

If the circulation is considered to be a periodic function then 

Pm - r 0 <*> ' i<Kt + 41 + 0U,) 

with 9 (z) being the phase angle of the circulation at some position, z, along 
the span, measured with respect to the value, ^) , at midspan of the plate so 
that 

I - pi e i(wt + ♦ 



r r H r 


H *1 


1 J 0 Pq(z)cos 0 dz + i J ( 


) l 0 (z) sin 8 dz I 



Assuming application of the lift force at the trailing edge the 
moment on the plate may be written as 
M = 2bN 

H 



M = 2 p Vbe 



P 



i(wt + § ) 



[/> 



(z)cos 8 dz + i I p^zjsin 9 dz I 

J o J ‘ 
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It is now assumed that 



P W 

\ 0 f 

where = circulation in the vortex 

u*j = the longitudinal fluctuating velocity component as measured by 
probe B 

and f = vortex shedding frequency. 

From reference (4) 

V* 

r 0 ~-f 

where V g = the velocity of the separation streamline. So that the assumption 
here is that 



u* ^ 



Substituting for p^(z) the moment may be written 



M rs/ pV be 



i(wt + 4 ) 



[a 



r H 

cos 9 dz + i 



1 1 0 4- sin0dz "J 



Now assuming the moment coefficient to be a real quantity and equating 
the two expressions for moment we have 



m^ rv 
o 



(fbHv] 



-1 



B 



cos 



9 dz 



1 



and 



r H — 



sin 9 dz = 0 



Here the assumption is that the wake is symmetrical behind the plate, so 
that the integral of an odd function is zero. Therefore, the forcing moment 
coefficient is proportional to the integral 



8 



£ 



u^ 2 cos 



6 dz 



where u^” is the longitudinal fluctuating velocity component measured by 
the traversing turbulence probe (probe B). 



B. Turbulence Probes 

Reference (4) describes the development and theory of the total head 
tube in detail. Briefly the following is a resume of the theoretical con- 
siderations involved in the measurements. 

In general the instantaneous response of a total head tube may be 
written 

u = i + El . (u + u') 2 (v + v* ) 2 (w + w') 2 

tg i y 2g 2g 2g 

where p is the pressure and u, v, and w are the velocity components in the 
coordinate direction. The barred quantities are mean values and the primed 
symbols are fluctuating values. 

From the geometry of the probe, it is a good assumption that the 
turbulence probe is much more sensitive to the velocity components parallel 
to the mean flow than to those which are perpendicular to the mean flow 
(i.e., the terms v ,2 /2g and w ,2 /2g are neglected). Also low frequency 
response of the piezoelectric crystal makes the turbulence probe insensitive 
to mean quantities. 

Therefore the instantaneous total head response becomes with these 
assumptions: 

u = £’ + 2uu» + uli 

tg % 2g 2g 

The D.C. portion of this signal (the time average of the fluctuation) 




u t2 

~2g 



is given by 
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The D.C* level is removed by a filter so that 

it2 



UU* 



U 1 



H tg " H tg + g + 2g 



uHE 

2g 



An important assumption at this stage is that all terms on the right 
hand side of the equation can be ignored with respect to uu*/g, whereupon 



H. - H. 
tg tg 



uuj 

g 



and the electrical signal as measured by the turbulence probe is proportional 
to the product of uu', the product of the mean longitudinal wake velocity and 
the longitudinal fluctuating velocity component. 



C. The Wave Analysis System 

The two electrical signals measured by the two turbulence probes located 
in the wake of the plate and positioned along the span were analyzed by the 
wave analysis system. 

In general the wave analysis system accepts one or two periodic or 
random signals as inputs and provides as output the power density of either 

i 

signal and the cross power density expressed as a complex number. The 
power density is a representation of the amount of power in the input signal 
contained in a given bandwidth. The cross power density is a measure of the 
power (contained in a frequency bandwidth) common to the two input signals. 

If the probe signals were harmonic they could be written as functions 
of time as f^(t) and fg(t) and their transforms as F^nj) and Fg(n 2 ). If 
F^(nj) is multiplied by Fg(n 2 ), the complex conjugate of Fg(n 2 ), and the 
product is averaged in time, only products of common frequencies remain so 

r (n) — — x (n) 

that i A g = F^(n)Fg(n) = cross power density. ^ ^ may also be 



written 
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§ AB (n) = \Tp iD A VPSDg 



This can be shovm since 



and 



therefore 






P,(n) - ^ J f,(t) e' inU ' t dt 



•h 
■h 



P 2 (n) - 5 






J f 2 (t) e' 1 ™’ 1 dt 






*, j r ’ 



h 



(t)f 2 (t +J-)dt - J f,(t)dt £ P 2 (n)e in “' (t +T) 
- t 4 ’ -h 



h » -c^ 



2j 

jf F^nje 1 " 40 ' 3 ' i j f,(t)e ln " , ' t 

at-co ^ 



dt 



-h 



1 

T, 



and 






J Z CO 

r,(t)r 2 (t +r)dt - \/^ n)F2(n, 3 « inu> ' J 






F,(n)F 2 (n) = = 



v\*» -w> 



A 



r 

J e in< 0 ,r dr^ J f,(t)f 2 (t + T )dt 

■H 



-h 



but 
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1 

T, 




e -ina>,? dT 



1 

T, 




f,(t)f 2 (t + 7 )dt 



and since machine operation gives 




F 2 (n) = \TpSD b 

then 




The TP 625 wave analyzer system accomplishes this by passing the random 
input signals f^(t) and fg(t) through a set of amplitude and phase matched 
band pass filters before multiplication. In order to preserve the phase 
angle 0 between the two input signals, two modes of operation of the 
multiplier yield the real and imaginary parts of a complex number called 
the co-spectrum and the quad-spectrum respectively so that 




co— spectrum) 2 + (quad- spectrum)^ 



^ -1 quad-spectrum 

co-spectrum 

For random signals these quantities have the units of power density or power 
divided by the filter band-pass at the center frequency of the filter. 

Now if it is assumed that the signal from probe B is correlated 
(i.e., similar) with the signal from probe A, the fixed probe, we can write 



and 

0 = 



U A = K(z) U B 



where K(z) represents the proportionality between the signals as a function 
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of their spanwise separation. As noticed previously the output of the 



probes are u’u so thatj^^J ^ u^u^u 2 . Note that it is tacitly 
assumed that u^ = Ug. Therefore 



14 J 



rv u'u'u 2 
A B 



U A U P 



and if we divide by 






l 2 we have 



\* 



AB 



u’u’u 2 * uju’u 2 
A B A B 






/v 
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AA 



u’uju 2 * u 2 
A A 



Substituting for uj* 



|4 



ABl 



A u 2 

Xaa 






K(z)u» s * K(z)u » 2 



K(z ) 2 u * 2 






Since wave analyzer operation gives 
^ = PSD. 



AA 



then 



14 J 



PSD. u 2 
A 






and referring to the discussion of the forcing moment coefficient where 
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cos 9 dz 



■ W 



and PSD, values 
A 



This suggests a method of evaluating the integral 
are available as outputs of the wave analysis system and values of u are 
easily measured in the tunnel. Values of 9 are available as shown from 
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^ -1 quad-spectrum 

co-spectrum 
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IV PROCEDURE 

A. Total Head Turbulence Probe Calibration 

The total head turbulence probes, electrometer preamplifiers, D.C. 
filters, and the turbulence probe signal amplifiers were calibrated as 
assembled units. The two system outputs were adjusted by changing the 
amplifier gain on one system so that the output of the two systems were 
equal for identical inputs. Figures 5 and 6, Appendix A, show the 
simplified turbulence probe signals when the probes were subjected to a 
step impulse. See reference (2) for calibration details. 

B. Tunnel Measurements 

The fixed turbulence probe was initially positioned in the tunnel at 
y/t = 0 and x/t = 1. See Figure 15, Appendix A for description of coordinate 
system. Alignment of the traversing turbulence probe (probe B) with respect 
to the fixed probe (probe A) is assured by the method of attachment of the 
„ two probes and the proper positioning of the traversing probe seal at the 
bottom of the test section. Measurements were taken with the plate free to 
vibrate at ambient velocities of 8, 11, and 14 fps, and clamped in a fixed 
position at an ambient velocity of 14 fps. The velocities were obtained 
from the contraction meter as described in reference (2). Mean wake 
velocities were measured as described in reference(3). While maintaining 
a constant ambient velocity by operation of the inlet and foot valves, the 
turbulence probe signals, plate vibration signal and a marker signal were 
recorded on four channels of the protable tape recorder for each spanwise 
position (z/t) of the traversing probe. Prior to and during the actual 
recording, the signals were monitored on oscilloscopes to insure that the 
proper signals were recorded. The pressure in the tunnel was maintained 
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at the same value by operation of the foot valve and was kept sufficiently 
high (3 psi) to prevent cavitation. 

C. Data Analysis 

The recorded signals were analyzed with the Technical Products Company 
TP 625 Wave Analyzer System. Each run consisted of a recording of the four 
signals (two fluctuating total heads, plate vibrational amplitude, marker) 
at a particular z/t and velocity. For each run the probe signals were 
rerecorded on a Sanborn-Ampex 2000 continuous loop tape recorder for con- 
tinuous play back into the wave analyzer system. Prior to each analyzation, 
the wave analyzer system was carefully calibrated with a 100 cps, O.46 volt 
peak to peak oscillator signal. The gains of the two wave analyzers were 
set the same and the oscillator (TP 626) output indicator was adjusted with 
a reference oscillator to ensure accurate frequency readings. The gains 
of the two analyzers TP 627 were matched and the Power Spectral Density set- 
tings of the multiplier TP 645 were adjusted to give the same output for 
„ identical inputs. In addition the phasing between the identical input 
signals was zeroed using the phase zero dial. Calibration in this manner 
ensures that the relative power levels in the signals are preserved in the 
output, and that zero phase angle is indicated for identical input signals. 
The output of the wave analysis system is in volts and is only a relative 
measure. A 50 cps. band width filter and a long averaging time (2 seconds) 
were chosen for the analyzation. The filter and 2 second averaging time 
were selected to give smooth curves from which data points could be 
accurately determined. 

For z/t = 1 and 9 at each velocity, the power spectrum density for 
each signal and the co-spectrum and quad-spectrum were obtained. For the 
remaining valves of z/t, the power spectrum density of M B W (the movable 
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turbulence probe), co-spectrum and quad-spectrum were obtained. 

The horizontal scale of the x-y plotter was adjusted so that a full 
sweep covered from 50 to 200 cps. The vertical scale of the x-y plotter 
was set to 1 volt per inch. Also for z/t = 9 at each velocity, the power 
spectrum density "B” was obtained using a 1 cps band width filter and an 
averaging time of 0.2 seconds to indicate the exact frequency of the 
energy in the wake. Similarly for the plate motion signal, a spectrum was 
obtained to pinpoint the frequency of plate vibration. 
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V PRESENTATION AND DISCUSSION OF RESULTS 

The outputs of the Wave Analyzer system are the power density spectrum 
of a single signal, the co-spectrum, and the quad-spectrum of the two signals 
taken jointly. These are presented for each run as Appendix C. 

For analyzation purposes, frequencies of 75 cps for runs with the test 
plate free to vibrate and 105 cps for runs with the test plate fixed were 
chosen. These frequencies were selected because they were the respective 
points of peak power as indicated by the plots of power spectrum density 
(PSD) and hence represent the frequency of vortex shedding. This reason 
for frequency selection is reinforced by the PSD of the vibrating plate 
signals which also has its peak power at 75 cps. See the power spectrum 
density for vibrating plate signals, runs 11, 16, and 21, Appendix C. 

Power density spectra were obtained from the wave analyzer system with 
narrow filter settings as well as with wide filter settings for runs 6, 11, 
16, and 21. These narrow filter power density spectra show more clearly 
the exact frequency of peak power than do the power density spectra made 
with wide filter settings. 

Table 1, Appendix B, lists the values of PSD, co-spectra and quad- 
spectra used in the following analysis. 

Power density spectra for probe A (fixed probe) are shown in Appendix C 
for runs 1, 6, 7, 11, 12, 16, 17, and 21. These are all runs with probe 
separation (z/t)= 1 and 9. It can be seen that for each velocity PSD^ is 
the same for all probe separations z/t. This indicates that at small probe 
separation, one probe does not affect the signal of the other probe and 
that the signal from the turbulence probe A is stationary. 

The power density for the fixed probe A increases in magnitude at the 
plate vibration frequency with an increase in ambient velocity. However 



18 



the power densities for the fixed plate (runs 1 through 6 of Appendix C, 
ambient velocity - 14 ft/sec) are higher than are the power densities for 
the vibrating plate runs with the same ambient velocity as the fixed plate 
runs (see runs 17 through 21 of Appendix C for vibrating plate runs with an 
ambient velocity of 14 ft/sec). Thus power density is not a function of 
ambient velocity alone. 

If one computes the self excitation parameter of reference (2), 

V* = — ^ j or run mac j e this study and compares and V*, 

it can be seen that decreases as the self excitation parameter V* 

increases. V* is the ratio of the transverse velocity of the trailing 
edge of the test plate to the ambient water velocity. 

It has been shown in Section III that the power density spectrum, 
is proportional to (u^u) 2 . Thus ^AA^ is proportional to (uj) 2 . The 
quantity, (u 1 ) 2 /^ may be considered a measure of the strength of the 
individual vorticies. Thus is proportional to the strength of 

vortices being shed. Figure 19 shows the strength of the vortices 
^^/u 2 f being shed plotted vs. the self excitation parameter V*. Figure 19 
shows that the vortex strength drops off rapidly with an increase in the 
self excitation parameter V*. Thus the increase in forcing moment co- 
efficient ra^ with an increase in the self excitation parameter V* cannot 
o 

be explained by an increase in the strength of the vortices being shed. 

For the case of the fixed plate (V* = 0) it appears that the values 
of vortex strength and coherency may vary with the ambient velocity. 

The power density spectra for the signal from probe B at small 
separations (runs 1, 7, 12, and 17 of Appendix C) are essentially the same 
as the power density spectra for the signals from probe A. 
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As shown in Section III: 

^ ^ Ab\ = ^( co ~ s P ec t rum ) 2 + (quad-spectrum) 2 ^ 

Q _ ^ aa ~l quad-spectrum 
co- spectrum 

Table 1 3.ists the values c^^ab! and 8 for a H runs. As explained in 
Section IV, Procedure, the wave analyzer system was adjusted so that the 
outputs were all relative to the 100 cycle sine wave calibration signal. 

The Technical Products Company wave analyzer system when in the power 
spectrum density mode of operation, multiplies the input signal by itself 
and in effect becomes a squarer. Time averaging then eliminates all 
products except those of common frequencies. Thus the power spectrum 
density of signal A = ^ . A ~ (u^u) 2 ; likewise ^ gg ^ (u^u) 2 . When the 
wave analyzer system is used to obtain the cross power spectrum density 
the two input signals are again multiplied and time averaged. 

When interpreting the machine output we must remember that all T s 
-are averages over a finite bandwidth. Thus when we look at ^5^ or ^>gg 
with the filter centered at frequency n Q , F^(n Q ) or Fg(n Q ) may be zero yet 
^ or can be non zero due to contributions at other frequencies 

within the filter pass-band. When looking at ^ A g we are most interested 
in its value at the vortex shedding frequency . If this frequency is only 
slightly different at the locations of probes A and B, and the random 
turbulent background at surrounding frequencies is negligibly small, 
can be zero and yet ^ and would be non-zero! 

Since |$ AB | ~ ( u a u b^ 2 )» <$AA ~ and ^ BB ~ (u B^ 2 ' then 

\$Ab\~ A A J ^$BB^ here U A and U B have the same f re 9 uencies ‘ 
Therefore^ is a measure of overall signal strength as well as a 

measure of the similarity of the two signal amplitudes. To eliminate the 
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effect of overall signal strength all values of |(J) for each velocity 
cam be normalized by their respective values (for each velocity ^ 

is stationary as shown in Appendix C). 

Therefore: ^ and is only a measure of the 

similarity of signal A and signal B. Values of are tabulated in 
Table 1, Appendix B. Figure 17 is a plot of vs. z/t for different 
values of the dimensionless self excitation parameter V* = 2b # 

Figure 18 is a plot of 0 vs. z/t for different values of V*. 

Vibrating plate nans at an ambient velocity of 1A fps were repeated 
after a two week interval to investigate the re produc ability of the 
experimental data. Figure 21 compares values of and 0 vs. z/t. It 
can be seen from this figure that the cross power spectra compare very 
favorably, but a phase angle variation is to be expected. 

From Figure 17 it can be seen that the correlation of the turbulence 
probe signals for vibrating plate runs decreases as probe separation 
increases which is what one would expect. Figure 17 also shows that the 
normalized cross power density increases for vibrating plate runs with an 
increase in the self excitation parameter V*. 

X ^ 

Figure 19 shows that the normalized cross power density (t) increases 
as the self excitation parameter V* increases. 

Figure 18 shows that for z/t< 5 for vibrating plate runs there is a 
gradual increase in © as z/t increases, but there is no significant change 
in 0 for an increase in the self excitation parameter V*. For probe 
separation z/t > 5 for vibrating plate runs 0 increases with an increase 
in z/t and decreases as the self excitation parameter V# increases. 

Figures 17, 18, and 19 show that the shedding of the vortices becomes 
more correlated and more in phase as the self excitation parameter V* 



increases. 
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For the self excitation parameter V* = 0 (fixed plate) with a velocity 

-r ■R' 

of 14 ft/sec. Figure 17 shows to decrease to near zero at z/t 7 and 
<3j> to increase for values of z/t > 7* For V* = 0 Figure 18 shows © to 
increase gradually for values of z/t 6 and 9 to increase rapidly to 
values y 180° for values of z/t y 6. This indicates that for V* = 0 
there may be a spanwise cellular structure with alternate cells 180° out 
of phase. Assuming symmetry the wavelength of the cell would be 14 body 
thicknesses. Roshko (6) reports a 18 diameter wavelength for measurements 
behind a circular cylinder at a body Reynolds number of 80. 

As shown in Section III: 



m ^ (bHfVu_ 2 ) 1 f — — — cos 0 dz 

0 §AA 

Values of cos 0 )/( ^ ^ ) were plotted vs. z/t and inte- 

grated from z/t a 1 to z/t = 9. It was assumed that the value of this 
integral is proportional to the integral evaluated over the entire plate 
span. 




was calculated by 




cos 0d(z/t) 



These values of m^ for the vibrating plate runs are plotted vs. V* in 

o 

Figure 20, Appendix 3. The dotted line in Figure 20 is a plot of m^. vs. 

o 

V* as determined from reference (2). It can be seen that these two plots 

of m^. are in qualitative agreement. In an effort to determine any pre- 

A3 | 2 / or Q in the integrand of the equation 

used to determine ra^ (i.e., is the amplitude coherence or the phase eo- 

o 

herence, the predominant factor in 



dominence of either 



1* 
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^ (fVu*) 1 



riy 

£aa 



COS 8 d(z/t)? ) it was assumed that cos 9 was 



constant, hence 



m ^ (fVu 2 )' 1 
r o 



/ 



s/t=9 



\M 

Iaa 



d(z/t) 



The resulting m^ plotted vs. V* is shown in Figure 20. 
o 



Next it was assumed that ( )/( ) was constant, hence 

-l r*A 

m^ /v» (fVu 2 ) j cos 9 d(z/t) 



m f ^ (fVu 2 ) 1 
o 



l 

The resulting m^ plotted vs. V* is shown in Figure 20. 
o 

The conclusion drawn from this fugure is that it is impossible to 
determine the predominate factor in 

- — cos 9 d(z/t) 

$ AA 

Undoubtedly both ( iU ! >/< ) and cos 9 are important. Indeed it may 

be said that amplitude correlation and phase correlation have no significance 

unless presented together. Phase correlation has no significance unless 

there is some amplitude correlation and vice versa. 

The • calculated values of m^. for the fixed plate runs assuming that 

o 

the integrals evaluated from z/t = 1 to z/t = 9 were proportional to the 
same integrals evaluated over the entire plate were not considered valid. 
However, by extrapolation of the integrand over the entire length, assuming 
a cellular structure with alternate cells 180° out of phase, values of 

o 

obtained are in qualitative agreement with reference (2). 

The co-spectrum of the vibrating plate runs, runs 7 - 21, Appendix C, 



show that some correlation exists at 140 cycles. This correlation was not 
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investigated, but it could be the correlation of the first harmonic of the 
vortex shedding frequency. The vibrating plate spectra for runs ((, ^ 2.1 
Appendix C, show power at about 140 cycles, and this could explain the 
correlation at 140 cycles. 
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VI SUMMARY AND CONCLUSIONS 



Spanwise cross correlations of the longitudinal turbulent velocity 
component in the early wake of a flat plate were measured in water under 
conditions of flow-induced torsional vibration of the plate in order to 
ascertain the nature of the known dependence of the forcing moment upon 
the ratio of vibratory plate velocity to free-stream fluid velocity. 

The following conclusions are drawn: 

1. The increase in the forcing moment coefficient m^ with the 

o 

self excitation parameter V* is due to an improvement in correlation 
or coherency of the wake along the span rather than to increasing 
strength of the individual vortices being shed. 



2. The forcing moment coefficient can be expressed in the form 






m f ^ (fbHVu 2 ) 
o 



cos 



9 dz 



0 $ AA 

3. Both the absolute magnitude of the cross power density 

and its phase angle 9 are important in explaining the relationship 

between the forcing moment coefficient m^ and the self excitation 

o 

parameter. Spanwise variations in phase angle are shown to be 
reduced with increasing V* and coherence of vortex strength at the 
vibrational frequency is increased. 



4. Measurements of the spanwise wake structure under conditions 
in which plate vibrations are prevented are consistent with a cellular 
structure with a phase shift of 180° between alternate cells. 

5. In order to minimize self-induced vibration, engineers should 
design structures in such a way as to eliminate wake coherency. This 
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might be accomplished by a random spanwise variation of trailing edge 
geometry, or by attaching to the trailing edge short independent 
sections which are free to rotate about the trailing edge. 
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VII RECOMMENDATIONS 

1. In the co-spectrum of probe B some indication of correlation was 
observed at 140 cps. It is recommended that this be investigated. 

2. It is recommended that measurements such as taken in this study 
be taken so as to include observations of the wake along the entire span. 
This will necessitate a new turbulence probe supporting mechanism. 

3. It is recommended that a 5, 10, or 20 cycle bandwidth filter be 
obtained for the Technical Products Wave Analyzer System. These more 
narrow filters will permit a more accurate determination of the frequencies 
of peak power in both the power spectrums and the cross power spectrums. 

The available one cycle filter proved to be too narrow for the random 
signal involved. 
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ILLUSTRATIONS 



TITL5S OF PICTURES ON THE FOLLOWING PAGE 



Figure 1: 
Figure 2: 

Figure 2 : 
Figure Z^: 



Water tunnel test section showing electrometers, 
amplifiers, DC filters, and strain gage controls. 



Left: Sanborn 2000 tape recorder. 

Center bottom: Precision Instrument Company portable 

tape recorder. 

Right: Technical Products wave analyzer system. 

Probes and supporting structure. 

Probes and supporting structure. 
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FIG. 5 

PROBE A SIGNAL 
CALIBRATION 



FIG. 6 

PROBE B SIGNAL 
CALIBRATION 




FIG. 7 

PROBE A SIGNAL 
VEL. = 0 FT. /SEC. 



FIG. 8 

PROBE B SIGNAL 
VEL. = 0 FT. /SEC. 




FIG. 9 

PROBE A SIGNAL 
VEL.= 11 FT. /SEC. 



FIG. 10 

PROBE B SIGNAL 
VEL.= 11 FT. /SEC. 
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FIG. 11 
PLATE SIGNAL 
VEL. = 0 FT. /SEC. 



FIG. 12 
PLATE SIGNAL 
VEL.= 8 FT. /SEC. 




FIG. 13 
PLATE SIGNAL 
VEL.= 11 FT. /SEC. 



FIG. 14 
PLATE SIGNAL 
VEL.* 14 FT. /SEC. 
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Itie origin of the coordinate system is located at the 
center of the test section at the trailing edge of the test 
plate . 



Figure 15: Turbulence probe arrangement, cross section 

of test plate, and coordinate system 
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FIG.16 BLOCK DIAGRAM OF INSTRUMENTATION FOR ONE or THE TWO 



TURBULENCE PROBES 
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Figure 17: Normalized cross power density vs. non-dimensional 

spanwise distance 
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? igure 18 : 



Correlation phase angle vs* non-dimensional spanwise 



distance 
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Figure 19: Normalized cross power density and vortex 

strength vs. self excitation parameter 
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Figure 20: Forcing moment coefficient vs. self excitation 

parameter 
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TABLJo 



Table 1 Measured Quantities 



Run 


z/t 

dimensionless 

distance 


PSD. 

A 

(volts)* 


co- 

spectrum 
(volts )* 


quad- 

spectrum 

(volts)* 


o< x 10 3 
(radians) 


u 

(ft/sec) 


V 

(ft/sec) 


1 


1 


3.25 


1.55 


0.12 


plate 


10.85 


14 


2 


3 


ft 


1.50 


0.20 


fixed 


ft 


ft 


3 


5 


If 


1.15 


0.33 


oc = o 


ff 


tt 


4 


7 


rt 


0.45 


-0.10 




tt 


ft 


5 


8 


it 


-0.55 


-0.17 




tt 


ff 


6 


9 


tt 


-1.25 


-0.73 




ft 


ff 


7 


1 


0.55 


0.31 


0.07 


3.93 


6.04 


8 


8 


3 


ft 


0.22 


0.06 


it 


tt 


tt 


9 


5 


If 


0.17 


0.05 


n 


ft 


tt 


10 


7 


ft 


0.10 


0.05 


n 


tt 


tt 


11 


9 


ft 


0.05 


0.05 


•t 


w 


tt 


12 


1 


1.63 


0.90 


0.15 


4.96 


9.50 


11 


13 


3 


it 


0.58 


0.30 


ft 


tt 


tt 


14 


5 


f! 


0.39 


0.25 


ft 


ff 


tt 


15 


7 


ff 


0.21 


0.20 


ft 


tt 


tt 


16 


9 


ft 


0.05 


0.20 


ff 


ft 


tt 


17 


1 


2.0 


0.83 


0.00 


4.63 


8.50 


14 


18 


3 


tt 


0.53 


0.21 


ff 


tt 


tt 


19 


5 


tt 


0.18 


0.24 


ff 


tt 


tt 


20 


7 


tt 


0.05 


0.17 


ff 


tt 


tt 


21 


9 


ft 


-0.10 


0.08 


ft 


tt 


tt 



* Values measured at 105 cps for fixed plate and 75 cps for vibrating plate 



Run 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13' 

14 

15 

16 

17 

18 

19 

20 

21 
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Table 2 Derived Quantities 



V* 

dimensionless 

velocity 


$ AB 
(volts) 


9 

(degrees) 


dimensionless 
power spectrum 


a IO 

u l 


1$ A6 \ cose 


UM/75 


plate 


1.55 


4.4 


0.48 


4.07 


6.260 


0.0198 


fixed 
V* = 0 


1.52 


7.5 


0.47 


3-98 


6.000 


ff 




1.19 


16.0 


0.37 


3.13 


3. 560 


ft 




0.46 


167.5 


0.14 


1.19 


-0.539 


ft 




0.58 


198.8 


0.18 


1.53 


-0.832 


ft 




1.45 


210.3 


0.45 


. 3. 81 


-4.770 


ft 


0.077 


0.32 


12.7 


0.58 


15.90 


4.950 


0.0150 


tt 


0.23 


15.3 


0.42 


11.50 


2.570 


ff 


it 


0.18 


16.4 


0.32 


8.80 


1.550 


ft 


tt 


0.11 


26.0 


0.20 


5.50 


O.546 


ft 


if 


0.07 


45.0 


0.13 


3.57 


0.173 


ff 


0.071 


0.91 


9.5 


0.56 


6.23 


5.560 


0.0181 


it 


0.65 


27.4 


0.40 


4.45 


2.570 


tt 


ti 


0.46 


32.7 


0.28 


3.11 


1.220 


u 


n 


0.29 


43.6 


0.18 


2.00 


0.417 


tt 


tt 


0.19 


73.2 


0.11 


1.44 


0.077 


tf 


0.052 


0.83 


3.6 


0.42 


5.83 


4.780 


0.0278 


it 


0.57 


21.7 


0.29 


4.02 


2.100 


if 


u 


0.30 


53.0 


0.15 


2.09 


0.375 


tf 


n 


0.18 


73.6 


0.09 


1.25 


0.063 


it 


it 


0.14 


146.3 


0.08 


1.11 


-0.113 


tf 



APPENDIX C 

POWER DENSITY SPECTRA 
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